Ab-initio simulation of photoinduced transformation of small rings in amorphous silica 
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We have studied the photoinduced transformation of small rings (3-membered) in amorphous silica 
by Car-Parrinello simulations. The process of ring opening leading to the formation of a couple 
of paramagnetic centers, namely an E' and a non-bridging-oxygen hole center (NBOHC), has been 
proposed experimentally to occur in silica exposed to F2 laser irradiation (at 7.9 eV). By using a 
new scheme for the simulation of rare events in ab-initio molecular dynamics (Iannuzzi, Laio and 
Parrinello, Phys. Rev. Lett. 90, 238303 (2003)), we have identified the transformation path for the 
opening of a 3-membered ring induced by a self-trapped triplet exciton, the migration of NBOHC 
and formation of a couple of stable E' and NBOHC paramagnetic defects. 

PACS numbers: 61.43.Fs; 61.72.Ji; 71.23.-k; 71.15.Pd 
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Extensive experimental and theoretical studies have 
been devoted to point defects in amorphous silica due 
to their relevance in the degradation of SiC>2-based elec- 
tronic devices and in the photosensitivity of optical fibers 
Q . For instance, it is well established that the possibility 
to change the refractive index and write Bragg gratings 
in optical fibers by UV illumination is connected to the 
photoinduced transformation of pre-existing defects or 
the formation of new ones. The formation of defects in 
amorphous silica upon F 2 excimer laser irradiation is also 
held responsible for the degradation of lenses for optical 
lithography in semiconductor technology |2(. 

In a recent experimental work Hosono et al. pro- 
posed that the main channel for color center formation 
in silica exposed to radiation of an F2 excimer laser (7.9 
eV, well below the Tauc band gap of a-Si02) is the gen- 
eration of an E' center and a non-bridging-oxygen hole 
center (NBOHC) from one-photon excitation and break- 
ing of strained bonds in small (3- or 4-membered) rings. 
The experimental evidence comes from combined Elec- 
tron Paramagnetic Resonance (EPR), Raman and optical 
spectroscopy measurements which show a correlation be- 
tween the intensities of the Raman lines associated to the 
breathing modes of 3- and 4-membered rings (D\ and D2 
lines) in samples with different fixation temperatures and 
the appearance of EPR and optical absorption signals as- 
signed to E' and NBOHC radicals. From these correla- 
tions Hosono et al. inferred that one-photon absorption 
processes at ~7.9 eV generate excitons self-trapped on 
the small rings which would lead to ring opening and 
pair formation of E' and NBOHC defects as shown in 
Fig. 1. 

In this paper we report the results of ab-initio simula- 
tions of the photoinduced processes leading to the open- 
ing of small rings in a-Si02 aiming at providing theoreti- 
cal support to the mechanism of E'-NBOHC pair genera- 




FIG. 1: Left panel: electron density of the hole of the 
triplet exciton self-trapped on a 3-membered ring in a 81- 
atoms model of a-Si02. The contour value on the electron 
density plot is 0.125 a.u.. The Kohn-Sham orbital is mostly 
localized on the oxygen atom 0(2). Right panel: path for the 
pair generation of the E' and NBOHC radicals upon photoex- 
citation of a 3-membered ring as proposed in Ref. 



tion proposed experimentally. The main outcome of the 
simulation is that the configuration of the two paramag- 
netic centers proposed experimentally and shown in Fig. 
1 is metastable on the excitate state, but the ring spon- 
taneously closes once the system is brought back on the 
electronic ground state. However, a stable E'-NBOHC 
pair is found via a further migration of the NBOHC. 

Ab-initio Car-Parrinello 0- d 01 molecular dynamics 
(MD) simulations have been performed within the frame- 
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work of density functional theory in the local spin den- 
sity approximation supplemented by generalized gradient 
corrections • Norm conserving pseudopotentials || , 
plane wave expansion of Kohn-Sham (KS) orbitals up to 
a kinetic energy cutoff of 70 Ry, a fictitious electronic 
mass of 800 a.u. and a time step of 0.15 fs have been 
used. The photoinduced transformations have been sim- 
ulated by adiabatic Born-Oppenheimer MD on the lowest 
triplet (Ti) excited state energy surface. Although the 
photoinduced reaction would occur on the singlet excited 
state, MD on the Ti excited energy surface is more eas- 
ily affordable within our framework and would provide 
transformation paths that also shed light on the pro- 
cesses induced by the singlet excitons as we have already 
demonstrated in a previous work on the photoinduced 
interconversion of oxygen deficient centers in a-SiC>2 0- 

Periodic models of a-SiC>2 containing 81 atoms at the 
experimental density of 2.2 g/cm 3 have been generated 
by quenching from the melt in classical MD, adopting the 
empirical pair potential by van Beest et ah [lQj . Quench- 
ing times as long as 5 ns have been used in the classi- 
cal simulations in order to generate models with a small 
number of 3- and 4-membered rings. Three models have 
been generated, one with one 3- and four 4-membered 
rings, a second with one 3- and five 4-membered rings, 
and a third one with no 3- and one 4-membered ring. The 
classical models have been then annealed at 600 K for 
1.5 ps by ab-initio MD. This procedure has been already 
demonstrated to provide models of a-SiC>2 with good 
structural, elastic and dielectric properties 
By exciting the system on the lowest triplet state, we 
have found that for the first among the models consid- 
ered the hole of the exciton is partially localized on the 3- 
membered ring. The self-trapped triplet exciton is shown 
in Fig. 1. Hereafter all the simulations presented would 
refer to this latter model. The triplet exciton induces a 
slight lengthening of the Si(l)-0(2) bond in Fig. 1 by 0.03 
A, but the ring is still locally stable. The transformation 
leading to the opening of the ring turns out to be an acti- 
vated process with an energy barrier much larger than the 
thermal energy. Thus it would not occur spontaneously 
on the time scale of Car-Parrinello simulations, few tens 
of ps long. In order to overcome this limitation, we have 
exploited a new technique recently devised to simu- 
late rare events within Car-Parrinello MD. The method 
is based on a coarse-grained dynamics in the space of 
few reaction coordinates, biased by a history-dependent 
potential which drives the system toward the lowest tran- 
sition state 0, 0, . The total energy of the locally 
stable configurations and transition states found along 
the transformation path have been then refined by con- 
strained geometry optimization 18]. The activation en- 
ergies of the different processes have been then estimated, 
although in principle also the entropic contributions to 
the activation free energies could be computed within the 
method of Ref. 0|. Due to the relatively large size of 



our simulation cell and the long simulation time needed 
to get accurate free energy estimates, we have restricted 
ourselves to the calculation of activation energies and not 
free energies. Entropic effects are small anyway for the 
breaking of siloxane bonds at room temperature as shown 
in Ref. . 

Following the scheme of Ref. [lj|j the collective re- 
action coordinates ^({R/}), function of the ionic posi- 
tions R/, define a set of associated collective variables 
s a which are treated as new dynamical variables. The 
extended system is described by the Lagrangian 

L = L ° +E^ M ^«-E^ fc 4^({ R i})- s -) 2 ( 1 ) 

a a 

- V(t,{s a }), 

where L a is the Car-Parrinello Lagrangian, the second 
term is the fictitious kinetic energy of the s a 's, the third 
term is a harmonic potential that restrains the value of 
the collective coordinates ^({Ri}) to the correspond- 
ing dynamic collective variables s a . V(t, {s a }) is the 
Gaussian-like history dependent potential defined in ref. 

0- 

In a first simulation, we have used a single collective 
variable which provides a measure of the length of the 
Si(l)-0(1) bond (cfr. Fig. 1) which is lengthened upon 
excitation on the triplet state. Namely, we have cho- 
sen the coordination number «5i(i)-o(2) between atoms 
Si(l) and 0(2). Following Ref. 14], the coordina- 
tion number between atom a and atoms b is defined as 
n a = J^b i ( r abyd where r a b is the distance between the 
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two atoms and scaling factor d is chosen as d=2.2 A. n a 
estimates the number of atoms b within the bond cut- 
off distance to atom a and decays smoothly for larger 
distances. Recent investigations have shown that the 
best efficiency of the method would be achieved in our 
case when the collective variable has a characteristic fre- 
quency {yfkjM) of the same order of magnitude as the 
Si-0 stretching mode 01- Accordingly, we have chosen 
M= 1.82 10 5 a.u. and k= 1 a.u. in Eq. (1) {y/k/M 
~ 10 14 s _1 ). The parameters which define the history- 
dependent Gaussian potential in Eq. (3) of Ref. 0] are 
5s x = 0.08 and W = 0.002 Hartree. MD simulations 
have been performed at 300 K [2oj |. 

Under the effect of the driving history-dependent po- 
tential the ring opens in 4.5 ps leading to the config- 
uration reported in Fig. 2b. Its energy, obtained by 
geometry optimization, is -0.27 eV with respect to the 
energy of the closed ring on the Ti excited state which is 
chosen hereafter as our zero of energy. At this point the 
electron of the triplet exciton is self-trapped on Si(l). 
In the configuration of Fig. 2b two paramagnetic de- 
fects are formed (an E' center on Si(l) and a NBOHC 
on 0(2), 3 A far apart) by breaking the Si(l)-0(2) silox- 
ane bond and overcoming an activation barrier of 0.38 
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FIG. 2: Snapshots of the locally stable configurations found 
along the path that leads to the formation of an E'-NBOHC 
pair from the photoinduced opening of a 3-membered ring. 
The atoms involved in the transformation are represented by 
spheres and connected by thicker sticks, (a) The 3-membered 
ring excited on the Ti state, (b) The adjacent E'-NBOHC 
defect pair, (c) The E' center and the 5-coordinated silicon 
atom, (d) The stable E'-NBOHC pair. 



eV. The energy barrier is obtained from constrained ab- 
initio MD simulations 0] at different values of a reaction 
coordinate chosen as the Si(l)-Si(2) distance. The con- 
figuration in Fig. 2b is locally stable on the Ti excited 
state, but by de-exciting the system on the ground state 
(S ), we have observed a charge transfer from the E' cen- 
ter to the NBOHC which leads to a spontaneous closure 
of the ring j^jj. Therefore, the two paramagnetic de- 
fects in the configuration of Fig. 2b are not stable in the 
ground state and cannot be assigned to the EPR signals 
observed experimentally. We may envisage to stabilize 
the E'-NBOHC pair by reducing the interaction between 
the two radicals to be realized by allowing a migration 
of one of two paramagnetic centers. In order to identify 
this process we have added a second collective variable in 
Eq. (1), defined as the coordination number rc-sifi)— 0(3) 
between atoms Si(l) and 0(3) in Fig. 2a. The presence 
of both the nsi(i)—0(2) an d n si(i)-o(3) collective vari- 
ables allows to explore different paths for the ring open- 
ing. The parameters M, k, 8s and W are the same as 
before. In this second simulation the system initially fol- 
lows the same path seen in the previous one which means 
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FIG. 3: Energies of the optimized local minima corrsponding 
to the configurations of Fig. 2. Closed ring, open ring, 5-fold 
silicon and E'-NBOHC pair refer to panels a, b, c and d of 
Fig. 2, respectively. The energy and geometry of the transi- 
tion states TS1, TS2 and TS3 are obtained from constrained 
geometry optimization [T3 | by using as reaction coordinate re- 
spectively the distance Si(l)-0(2), Si(l)-0(3) or Si(4)-NBO 
in Fig. 2. 



that the Si(l)-0(2) bond is weaker than the Si(l)-0(3) 
one. However, as the simulation proceeds, the history- 
dependent potential drives the system away from the lo- 
cal miminum of Fig. 2b. First, the Si(l)-0(2) bond is 
recovered and then a new conformation of the open ring 
is formed by breaking the Si(l)-0(3) bond. The NBOHC 
generated in this way (on 0(3)) floats around and finally 
forms a new bond with Si (4) which becomes 5-fold co- 
ordinated. This process has a higher activation barrier 
(0.64 eV) but leads to a configuration (Fig. 2c) which is 
locally stable with a total energy of -0.54 eV. From this 
configuration it is now envisageable to generate a new 
NBOHC by breaking one of the five Si(4)-0 bonds. To 
identify which Si(4)-0 bond is most prone to break, we 
have started a third simulation from the configuration of 
Fig. 2c with the total coordination number nsiU)—o as 
collective variable. The final structure of this latter pro- 
cess is shown in Fig. 2d. Its total energy is -0.28 eV and 
the activation barrier for the formation of the NBOHC is 
0.46 eV (Fig. 3). The E' and the NBOHC centers in this 
final configuration are distant enough (6.3 A) to prevent 
the charge transfer upon de-excitation on the S c state. 
As the two radicals do not interact each other, there is 
no singlet-triplet splitting and the EPR signals are sup- 
posed to be equal to those of the isolated paramagnetic 
centers. The path identified by the local minima of Figs. 
2a-c represents thus a viable mechanism for the open- 
ing of the ring and migration of the NBOHC suitable to 
provide a pair of stable paramagnetic defects. The total 
trajectory starting from the initial configuration of Fig. 
2a and reaching the final configuration of the stable E'- 
NBOHC pair in Fig. 2d is 2.7 ps long . The energies of 
the configurations in Fig. 1 and of the transition states 
separating the different local minima are collected in the 
scheme of Fig. 3. The geometry of the transition state 
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FIG. 4: Geometry of the transition state TS2 of Fig. 3. The 
distances between the silicon atom Si(l) and the oxygen atoms 
0(2) and 0(3) are 1.80 A and 2.05 A, respectively. The tetra- 
hedron formed by atom Si(l) is highly distorted and under- 
goes a puckering once a NBOHC is formed on 0(3) (cfr. Fig. 
2c). As a consequence the resulting E' center that does not 
point towards the NBOHC. 



leading to the 5-fold coordinated Si is reported in Fig. 4. 
We remark that the recovery of the Si(l)-0(2) bond mov- 
ing from configuration in Fig. 2b to the transition state 
in Fig. 4 does not restore the initial configuration of the 
closed ring (Fig. 2a) . No intermediate local minima from 
Fig. 2b to Fig. 2c is found. The activation barrier for the 
whole process (the energy of the highest transition state) 
is 0.64 eV, sufficiently low to make this transformation 
channel viable at room temperature. These values for the 
activation energies are supposed to depend on the local 
environment and strain of the small ring. As a conse- 
quence only a fraction of the small rings are expected to 
undergo this trasformation upon photoexcitation. 

In summary, we have identified the transformation 
path for the opening of a 3-membered ring induced by 
a self-trapped triplet exciton. Once open in the ex- 
cited electronic state, the 3-membered ring gives rise to 
a couple of strongly interacting E' and NBHOC defects, 
which recombine quickly after electronic de-excitation. 
However, a simple path with a low energy barrier (0.64 
eV) has been identified for the migration of the NBOHC 
which leads to the formation of a couple of stable E' 
and NBOHC paramagnetic centers. The simulations thus 
provide a theoretical support to the interpretation of the 
experimental data in Ref. [2| on the formation of color 
centers in a-Si02 by F2 laser irradiation. 
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